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Generalized parton distributions (GPDs) have become a standard QCD tool for analyzing and 
parametrizing the non perturbative parton structure of hadron targets. GPDs might be viewed as 
non-diagonal overlaps of light-cone wave functions and offer the opportunity to study the partonic 
content of the nucleon from a new perspective, allowing one to study the interplay between longi- 
tudinal and transverse partonic degrees of freedom. In particular, we will review some of the new 
information encoded in the GPDs through the definition of impact-parameter dependent parton 
distributions and form factors of the energy-momentum tensor, by exploiting different dynamical 
models for the nucleon state. 
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1. Introduction 

Generalized Parton Distributions (GPDs) have become a standard QCD tool for analyzing and 
parametrizing the non perturbative parton structure of hadron targets, for reviews see |[T], ^, ||, ^ |5|]. 
GPDs have been introduced in the past in different contexts (see, e.g., |Q, ^), but have raised a large 
interest in the hadron community only when their importance was stressed in studies of deeply 
virtual Compton scattering (DVCS) ^ ^ |l^] and hard meson production [ [Til ] in connection with 
the possibility of factorizing their contribution and gaining information on the spin structure of 
the nucleon Being defined in terms of nondiagonal matrix elements of the same correlation 
functions entering the definition of the parton distribution (PDs), GPDs reveal the partonic content 
of hadrons from a complementary perspective. They do not represent any longer a probability, but 
rather the interference between ampUtudes describing different parton configurations of the nucleon 
so that they give access to momentum correlations of partons in the nucleon. Furthermore, the finite 
momentum transfer to the proton makes a second space-time structure of the process possible, and 
after the Fourier transform in the impact-parameter space allows one to define spin-dependent 
densities which describe how partons are spatially distributed in the transverse plane [p^]. 

The GPDs can also be viewed as the generating functions for the form factors of the twist- 
two operators governing the interaction mechanisms of hard processes in the deep inelastic regime. 
These generalized form factors do not couple directly to any known fundamental interactions, but 
can be studied indirectly looking at moments of the GPDs. The most peculiar example are the form 
factors of the energy-momentum tensor, which give information about the spatial distribution of 
energy, angular momentum and forces experienced by quarks and gluons inside hadrons. 

After summarizing in sect. 2 definition and basic properties of GPDs, in sect. 3 we show 
results for the GPDs in the impact-parameter space using a light-cone quark model and discussing 
in particular the correlations of spin and orbital angular momentum of the quarks in the nucleon. 
Finally, in sect. 4 we review some results for the form factors of the energy-momentum tensor in 
the framework of the chiral-quark soliton model. 

2. Definition and basic properties of GPDs 

Parton distributions are defined in terms of matrix elements of light-cone bilocal operators 
between proton states of equal momenta. In general, with initial (final) momentum p {p') and 
helicity Aat (A^) one defines a set of quark generalized quark distributions for a hadron with spin ^ 

(y,A;|€?^(x,0^)|/j,A;v), (2.1) 

with 

^^(^,0^) = J ^e'^'^^-Y,{-'-,0^)rY/-,0^). (2.2) 

In Eq. ( |2.2[ ) P = {p + p')/2, and the operator F is a matrix in Dirac space which selects different 
spin polarizations of the quark fields. For three particular matrices F one can classify eight lead- 
ing twist GPDs: /) two unpolarized quark GPDs, H'^ and E'', for F = 7+; //) two longitudinally 
polarized quark GPDs, //^ and E'^', for F = 7^75; Hi) four quark chiral-odd GPDs for F = ia'^js, 
i.e. Hj- and Ej-, which involve the density operator for transversely polarized quarks, and Hj- and 
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Ej, defined in terms of a quark operator whicli fiips the transverse spin of tlie quark. Tiie tilded 
distributions correspond to matrix elements between nucleon states with flip of the polarization, 
while the untilded distributions refer to no-flip of the nucleon polarization. Analogous definitions 
hold for the gluon GPDs. 

Because of Lorentz invariance the eight GPDs can only depend on three kinematical variables, 
i.e. the (average) quark longitudinal momentum fraction x = /P^ , the invariant momentum 
square ? = = {p' — p)^, and the skewness parameter ^ given by i§ = — A+/(2P+). In addition, 
there is an implicit scale dependence in the definition of GPDs corresponding to the factorization 
scale /i^ used to separate the (universal) matrix element defining a GPD inside the entire amplitude 
describing the process under study. 

In the forward case, p = p' , both A and t, are zero. In this case the functions W^, W and Hj 
reduce to the usual DIS parton distribution functions, i.e. the quark density, helicity and transversity 
distributions, respectively. No corresponding relations exist for the functions E'^, Ej. and Hj-, 
because in the forward limit they decouple in their defining equations. However, they do not vanish. 
In particular, E''{x,0,0) carries important information about the quark orbital angular momentum. 
In contrast, £'|(;c,0,0) vanishes identically being an odd function of by time reversal symmetry. 

Moments in the momentum fraction x play an important role in the theory of GPDs. Weighting 
Eq. (lA) with integer powers of x and integrating over x, the correlation function ff^ reduces to a 



local operator and the corresponding matrix elements can be parametrized in terms of generalized 
form factors (GFFs) related to Mellin moments of the GPDs. More specifically, one has 
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for the unpolarized GPDs; 
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for the polarized GPDs; 
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odd 



for the chiral-odd GPDs. 

For the lowest moment « = 1 in Eqs. (U)-(||4l) one finds A\ ^{1) = Ff (t) , B%{t) = F^{t), 
^1 o(0 =§1(0 andBj ()(?) = gp{t) where Fj^, and gp are the quark contribution to the Dirac, 

Pauli, axial and induced pseudoscalar form factors, respectively. Furthermore, in the chiral-odd 



3 



Generalized Parton Distributions 



B. Pasquini 



sector one finds A|.j ^{t) = gj{t) and 2A^j q(0) o(0) = where gj is the quark tensor form 
factor and K"| describes how far and in which direction the average position of quarks with spin in 
the x-direction is shifted in the j?-direction in an unpolarized nucleon [[l3|]. 

The second Mellin moments of unpolarized GPDs can be related to the form factors of the 
energy-momentum tensor (EMT) of QCD by [ |T4| ] 

A%{t) = Ml{t) + \dl{t)E,\ 4o(0 = 2J^{t) -Ml{t) - \d'i{t)£,\ (2.7) 

The form factor M2{t) at f = reduces to the second Mellin moment of unpolarized parton distri- 
butions accessible in inclusive deep inelastic scattering, and represents the fraction of the nucleon 
momentum carried by quarks. The form factor d\{t) provides information on the distribution of 
strong forces in the nucleon, similarly as the electromagnetic form factors contain information 
about the electric charge distribution. The form factor J'^{t) is relevant for the spin structure of the 



nucleon thanks to the so called Ji's sum ruleHa, 10] 



J'i{t = 0) = {f\=S' AlJO)+B%{0) 



(2.8) 



where {J'^) is the total angular momentum along the direction / carried by quarks and antiquarks in 
a proton with spin S'. In the case of a proton polarized in the positive £-direction, one can further 

split the Ji's sum rules into spin and orbital angular momentum parts, i.e. - = -Z^ + where the 
contribution from the quark spin can be obtained from the moments of the usual polarized quark 
densities, i.e. If^ = A\ ^{t = Q) . The orbital angular momentum has recently been calculated in 
lattice simulations by the LHPC [jT5|] and QCDSF [16] collaborations. These calculations are for 



pion masses as low as 350 MeV and volume as large as (3.5 fm)^, providing results in the MS' 
renormaUzation scheme. Extrapolation to the physical pion mass requires a combination of full 
QCD lattice and Chiral Perturbation Theory [p^]. Two remarkable features are found. The first 
is that the magnitude of the orbital angular momentum contributions of the up and down quarks 
are sepai^ately sizable, L" f» — L'' k, 0.30, yet they cancel nearly completely at all pion masses, 
jji+d ~ [15, |l8|], indicating that the total angular momentum of quarks in the nucleon is of the 
same size as the quark spin contribution. The second is the close cancellation between the orbital 
and spin contributions of the down quarks for all pion masses, J'' ^0 [ [Ist |T8| ]. However, before 
drawing definite conclusions, one should be aware that these results do not include the contributions 
from disconnected graphs. Such contributions cancel in the difference of u and d quark distributions 
but may well be important in their sum. 

From the experimental side, first model-dependent constraints on the angular momentum were 



extracted from recent DVCS data. The analysis was performed comparing data from HERMES 019| ] 



on transverse-target and beam-charge asymmetries, and cross section data from JLab [gOp with 



various GPD-models having as free parameters [^, |22|]. Although the extracted values are 
strongly model-dependent, such analysis show for the first time that DVCS data have indeed the 
potential to provide quantitative information about the spin content of the nucleon. 



Whereas Eq. ( |2.8| ) provides the angular momentum carried by the quarks regardless of their 
spin, one can also investigate how much each quark polarization component contributes to {J'^). 
Taking for example the case of transverse polarization in the x direction, one can decompose 
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with respect to quarks of definite transversity, i.e. = -/^ + Jq _f where 7^ _|_. corresponds to 
the angular momentum in the x direction carried by quarks with transverse polarization in the iti- 
direction The transversity components Jq^^ can be related to the second Mellin moments of 
GPDs as follows 



5^ 



(2.9) 

In the case of transversely polarized nucleon, -Z^+f + gives the sum rule (^!^), while for 
unpolarized target only the second term in Eq. contributes. Thanks to this decomposition 
of the angular momentum in transversity components one can gain information on the correlation 
between the transverse spin and the transverse angular momentum carried by the quarks in an 
unpolarized target. Note that the same linear combination of GPDs {2Hj{x,0,0) +£'7-(;c,0,0)) that 
appears in Eq. ( |2.9| ) also describes the transverse displacement of quarks with a given transversity 
in an unpolarized target relative to the center of momentum (see sect. 3). 

3. Spin densities in the impact-parameter space 

A convenient way to make explicit which kind of information on hadron structure is contained 
in the GPDs is the representation in terms of overlaps of light-cone wave functions (LCWFs) which 
are the probability amplitudes to find a given A^-parton configuration in the Fock-space expansion 
of the hadron state. In the following, we will confine our analysis to the three-quark sector, by 
truncating the light-cone expansion of the nucleon state to the minimal Fock-space configuration 
and adopting a light-cone constituent quark model (CQM) which has been successfully applied in 



the calculation of the electroweak properties of the nucleon []23|]. 

As outlined in Ref. [^], the starting point is the three-quark wave function obtained as solution of 
the Schrodinger-like eigenvalue equation in the instant-form dynamics. The corresponding solu- 
tion in light-cone dynamics is obtained through the unitary transformation represented by product 
of Melosh rotations acting on the spin of the individual quarks. In particular, the instant-form wave 
function is constructed as a product of a momentum wave function which is spherically symmetric 
and invariant under permutations, and a spin-isospin wave function which is uniquely determined 
by SU(6) symmetry requirements. By applying the Melosh rotations, the Pauli spinors of the quarks 
in the nucleon rest frame are converted to the light-front spinors. The relativistic spin effects are 
evident in the presence of spin-flip terms in the Melosh rotations which generate non-zero orbital 
angular momentum components and non-trivial correlations between spin and transverse momen- 
tum of the quarks. On the other hand, the momentum-dependent wave function keeps the original 
functional form, with instant-form coordinates rewritten in terms of light-front coordinates. Model 
results for the non-polarized, polarized and chiral-odd GPDs in the momentum space have been 



presented in Ref. [ |24| , [25| ], with an extension of the model to include the contribution from higher- 



Fock state components in Ref. [|25|]. Here we focus on the GPDs in the impact-parameter space, 
which give complementary information to the transverse momentum dependent parton distributions 
on the spin-spin and spin-orbit correlations of quarks in the nucleon p^]. 

When ^ = and x > 0, by a two-dimensional Fourier transform to impact-parameter space, 
GPDs can be interpreted as densities of quarks with longitudinal momentum fraction x and trans- 
verse location b with respect to the nucleon center of momentum [|^, Depending on the 
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Figure 1: The spin-densities for (transversely) i-polarized quarks in an unpolarized proton (left panels) and 
for unpolarized quarks in a (transversely) i-polarized proton. The upper (lower) row corresponds to the 
results for up (down) quai^ks. 



polarization of both the active quark and the parent nucleon, according to Refs. [27, 28] one de 



fines tiiree-dimensional densities representing the probability to find a quark with transverse spin 
St in the nucleon with transverse spin Sj- It reads 

1 



p{x,b,ST,ST) = \ 



H{xy)+s'TS'T{HT{xy) 



4M2 



AhHrixy 



+^ (s'TE'{x,b^)+s'T \E'T{xy)+2H!rix,b^)]) U^^yy -b^5ij)s{^H'i{x,b^) 



,(3.1) 



f. In Eq. (3.1) enters the 



M 

where the derivatives are defined /' = ^ /, and A^f = 4-^ 
Fourier transform of the GPDs in the impact-parameter space. 
In Eq. (p|), apart from the two orbitally symmetric monopole terms in the first line, there are 
two dipole structures proportional to b^e^'s'j. and b^e^'Sj, and a quadrupole term proportional to 
s'j-{2b'bj -b^6ij)s{. The (derivatives of the) GPDs E{x,b^), ET{x,b^) +2HT{x,b^) and HT{x,b^) 
thus determine how strongly the orbital symmetry in the transverse plane is distorted by the dipole 
and quadrupole terms. 

Lattice calculations accessing the lowest two x-moments of the quark transverse-spin densities 
have recently been presented in Ref. [29|. Here we show some results in the light-cone CQM 
for the first x-moment of the spin distributions in the cases of transversely polarized quarks in an 



unpolarized nucleon and unpolarized quarks in a transversely polarized nucleon, referring to ||30p 
for the discussions of more complex spin configurations. 

In the case of transversely polarized quarks in an unpolarized proton the dipole contribution 
Ej{x,b^) +2Hj{x,b^) introduces a large distortion perpendicular to both the quark spin and the 
momentum of the proton, as shown in the left column of Fig. [T[ Evidently, quarks in this situation 



also have a transverse component of orbital angular momentum. This effect has been related [13] 
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to a non-vanishing Boer-Mulders function which describes the correlation between intrinsic 
transverse momentum and transverse spin of quarks. Such a distortion reflects the large value of 
the anomalous tensor magnetic moment Kj for both flavors. Here, Kj = 3.98 and = 2.60, to be 



compared with the values Kj and K^j 1.9 of Ref. [29|. Since Kj ~ —h^, the present results 
confirm the conjecture that is large and negative both for up and down quarks [|l3|]. 

As also noticed in Refs. [|l^, the large anomalous magnetic moments k:" '^' are responsible 
for the dipole distortion produced in the case of unpolarized quarks in transversely polarized nucle- 
ons (right column of Fig. [l]). With the present model, k" = 1.86 and = —1.57, to be compared 
with the values k" = 1.673 and k'^ = —2.033 derived from data. This effect can serve as a dynami- 
cal explanation of a non-vanishing Sivers function which measures the correlation between the 
intrinsic quark transverse momentum and the transverse nucleon spin. The present results, with the 
opposite shift of up and down quark spin distributions, imply an opposite sign of f^j for up and 



down quarks ||31|] as confirmed by the recent observation of the HERMES [g2|] and COMPASS [|33|] 



collaborations. The results in Fig. H are also in qualitative agreement with lattice calculations [^]. 

4. Energy-Momentum Tensor 

The nucleon EMT form factors were studied in lattice QCD [[I^], chiral perturbation the- 



ory [|34|], and models such as the Skyrme model [ ]35| ] and chiral quark soliton model (;^QSM) [ ]36| , 



37[ ]. Here we summarize some interesting features investigated within the ;|^QSM when looking 



at the spatial distribution of the EMT form factors [|l^, 36]. The ;^QSM provides a field theoretic 
description of the nucleon in the limit of a large number of colors Nc, where the nucleon appears as 
chiral soliton of a static background field. Numerous nucleonic properties have been described in 
this model, giving predictions in agreement with phenomenology within an accuracy of (10-30)%. 
The spatial distribution of the EMT form factors can be obtained from the Fourier transform with 
respect to A of the matrix element of the static EMT calculated in the Breit frame [|l^]. The nor- 
malized energy density Aiir^pE{r) /M^ is shown in Fig. §a as a function of r in the chiral limit of 
a vanishing pion mass and for the physical situation with a pion mass of 140 MeV. In the center 




Figure 2: (a) The normalized energy density 4-Kr^pE{r)/MN from the ;|^QSM as a function of r in the chiral 
limit of niji = (dashed curve) and for nij^ = 140 MeV (solid curve), (b) The same for the normalized angular 
momentum density 47i:r^pj(r)/7Ar (taken from |36|). 
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Figure 3: (a) The pressure p{r) from the ;i^QSM as function of r for = (dashed curve) and 140MeV 
(solid curve), (b) r^p{r) as function of r at the physical value of (taken from [^). 



of the nucleon the energy density is Pe{0) = 1.70 GeV/fm^ or 3.0 x lO^^ g cm ^, corresponding 
roughly to 13 times the equilibrium density of nuclear matter. As the pion mass decreases, the 
energy density is spread more widely. According to the role of the pion field in the ;^QSM, where 
one can associate the contribution of the discrete level to the quark core and the contribution of the 
negative continuum states to a pion cloud, this means that the range of the pion cloud increases and 
the nucleon becomes larger. Actually, the mean square radius increases from 0.67 fm^ in the 
case of the physical pion to 0.79 fm^ in the chiral Umit. With increasing pion mass up to 1.2 GeV 
this trend is confirmed with the nucleon becoming smaller and smaller [06|]. 

The normalized radial distribution of angular momentum 4Kr^pj{r)/JN, is shown in Fig. ^ 
as a function of r for mji = and 140MeV. For any nij^ at small r one finds pj{r) r^. The mean 



square radius [rj) decreases with increasing niji [36] in agreement with the idea of a shrinking pion 
cloud. For a physical pion one finds {rj) = 1.32 fm^. At large r in the chiral limit pj{r) l/r^ 
such that (rj) diverges. 

Fig. |3|a shows the pressure p{r) as function of r. In the physical situation p{r) takes its global 
maximum at r = with p{0) = 0.23GeV/fm^ = 3.7 • 10^"* Pa. This is i^(lO-lOO) higher than 
the pressure inside a neutron star. Then p{r) decreases monotonically (becoming zero at ro = 
0.57 fm) till reaching its global minimum at rp = 0.72 fm, after which it increases monotonically 
remaining, however, always negative. The positive sign of the pressure for r < ro corresponds to 
the repulsion among quarks imposed by Pauli principle, while the negative sign in the region r > ro 
means attraction in agreement with the idea of a pion cloud responsible for binding the quarks to 
form the nucleon. The subtle balance between repulsion and attraction, ultimately producing a 
stable soli ton, can be better appreciated from Fig. ^ showing r^p{r), where the shaded regions 
have the same surface areas but opposite sign and cancel each other within numerical accuracy. 
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